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Nomenclature 

COM 

center of momentum 

da 

7TE 

[Omar j3 

spectral distribution, 2?r pi dO F. ((sin 0 

J 0 dp 3 

E 

pion total energy 

Ecm 

center of momentum energy 

e £± 

dp 3 

LIDCS 

EAS 

extensive air shower 

LIDCS 

Loren tz invariant differential cross section 

m p 

proton mass 

777 7T 

pion mass 

Pp 

proton momentum 

V 

pion momentum 

Pmax 

maximum possible momentum scattered pion can have for given sfs 

P± 

pion transverse momentum, psin 0 


magnitude of center of momentum frame four momentum, equal to total energy in 
center of momentum frame 

T 

pion kinetic energy 

^lal) 

laboratory frame kinetic energy of incoming proton 

0 

angle of pion scattering with respect to direction of incident particle 

u 

max fPn lax ,/3 _ r 2 n f) 

total cross section, 2tt / dO I dp E / * 

^ ^ Pmm dp \Jp* - h 

Quantities with an asterisk (e.g., 0*) refer to the quantities in the center of momentum frame 
whereas quantities without an asterisk (e.g., 0) refer to the quantities in the laboratory frame. 
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Abstract 


An accurate knowledge of cross sections for pion production in 
proton-pwion collisions finds wide application in particle physics , 
astrophysics, cosmic ray physics, and space radiation problems, 
especially in situations where an incident proton is transported 
through some medium and knowledge of the output particle spectrum 
is required when give n the input spectrum. In these cases, accurate 
parameterization# of the cross sections are desired. In this paper 
much of the experimental data arc reviewed and compared with a 
wide variety of different cross section parameierizations. Therefore , 
parameierizations of neutral and charged pion cross sections are pro- 
vided that give a very accurate description of the experimental data . 

Lorentz invariant differential cross sections, spectral distributions, 
and total cross section parameierizations arc presented , 

1. Introduction 

Pion production in proton-proton collisions has been extensively studied for many years; this 
knowledge now finds useful applications in a variety of areas as follows: 

1. Two important types of particle detectors are the hadronic and electromagnetic calorime- 
ters (ref. 1) where an electromagnetic or hadronic shower is initiated by a liigli -energy 
incoming particle; from a Monte-Carlo simulation of the shower, one is able to deduce 
important characteristics of the incoming particle such as its energy and identity 

2. The primary cosmic rays can be detected by a variety of methods, depending on the 
incident energy; for the very high-energy cosmic rays, where the flux is relatively low, the 
extensive air showers (EAS's) (refs. 2, 3, and 4) provide the most convenient means of 
detection; the EAS is analogous to the hadronic or electromagnetic calorimeter used in 
particle physics but with the atmosphere of the Earth being the active volume in which 
the shower develops; the EAS has both ..electromagnetic and hadronic components, and 
similar to the calorimeter, the energy and identity of primary cosmic ray nuclei can be 
deduced via. Monte-Carlo simulation of the showers (refs. 2 and 3) 

3. In long-duration human space flights, such as a mission to Mars, the radiation levels 
induced by galactic cosmic rays can exceed exposure limits set for astronauts (refs. 5 
and 6); in determining the radiation environment inside a spacecraft, one needs to 
transport the exterior cosmic ray spectrum through the spacecraft wall to determine the 
interior radiation spectrum 

4. In gamma ray (refs. 7 and 8) and high-energy neutrino astronomy (refs. 9 and 10), 
the diffuse background radiation is due in large part to the gamma rays and neutrinos 
produced in proton collisions with the protons in the interstellar medium; in addition, pion 
product ion from proton-proton collisions finds applicat ions in the calculation of gamma 
ray emission from the accretion disk around a black hole (ref. 11) 

In all these applications, having an accurate knowledge of the cross sections for pion 
production in proton-proton collisions is crucial. In addition, most of the applications mentioned 
require solving the transport equations that determine the particle spectrum on one side of a. 
material (active volume of calorimeter, atmosphere of Earth, spacecraft wall, or interstellar 
medium) given the incident part icle spectrum. Use of pion production cross sections in such 
transport codes requires that the cross section be written in a simple form. The transport codes 
have many iterative loops, which will take too much computer time if the cross section formulas 
also contain many iterative loops. Thus it is most advantageous if one can write down simple 



formulas which parameterize all experimental data on pion production cross sections; this is the 
aim of the present work. 

In this paper, simple algebraic parameterizations of charged and neutral pion production 
cross sections valid over a range of energies are presented. The cross sections provided are 
Lorentz invariant differential cross sections (LIDCS’s), laboratory frame spectral distributions 
(i.e., energy differential cross sections), and total cross sections because they are the types of 
cross sections most widely used in transport equations. Many such parameterizations have 
been presented before, but the problem is deciding which are correct and whether a particular 
parameterization applies only to a limited data set or is valid over a wider range. In the present 
work, an exhaustive data search has been performed, and as many different parameterizations as 
possible have been compared with as much data as possible so that definitive conclusions could 
be reached concerning which is the most accurate parameterization to use. 

The cross sections discussed in this p>aper are for inclusive pion production in proton-proton 
collisions; that is, the reactions considered are p +p — * tt+ A', where p represents a proton, 
7T represents a. pion, and X represents any combination of particles. An extensive search for 
LIDCS data was performed, and the data were used to compare all available parameterizations. 
A method for generating parameterizations for these cross sections is also described and applied 
to 7T^ production. Spectral distribution and total cross section formulas were not developed 
directly because of lack of data. Instead, the most successful LIDCS parameterizations were 
first transformed into laboratory frame spectral distributions by numerical integration. These 
spectral distributions were parameterized and then numerically integrated to generate laboratory 
frame total cross sect ions. Finally, the total cross sections w 7 ere compared with available data and 
parameterized as well. This procedure is discussed, and (lie parameterizations of the numerical 
results are given. Multiple checks of the accuracy of all results were made, and some of them 
are presented. 

Finally, there are a lot of figures in this paper. What often happens when various authors 
come up with a parameterization is that they only apply it to a limited data set. Often when a 
particular parameterization is applied to other data, it does not work; this is the reason for the 
large number of figures in the present paper. The aim is to show 7 that the parameterizations in 
this paper do apply w r ell to a whole range of experimental data. 

The derivation of the maximum momentum is presented in appendix A. The kinematic 
relations between the two reference frames- the center of mass and the laboratory frames — 
are presented in appendix B. Appendix C presents a synopsis of data transformat ioas. 


2. Comparison of Lorentz Invariant Differential Cross Sections 


The object is to determine an accurate parameterization for inclusive LIDCS's, wdiich can 
be confidently applied to regions where no experimental data are available. For example, if 
the formulas were to be used for the purpose of developing radiation shielding materials, the 
parametric equation would need to be extrapolated to energies lower than those for wdiich data 
are available. The most convenient formulas are those that are in closed form, since they are 
easily used and take relatively little computer time in numerical calculations. Some of the 
formulas that were considered as representations of the LIDCS s were not in closed form but 
included tabulated functions of energy (i.e., numerical values were given for specific energy values 
rather than a functional form). When comparing parameterizations, closed-form expressions 
were given precedence over other equally accurate formulas. 


The invariant single-particle distribution is defined by 
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where 


£sl 

dpi 


is the differential cross section 


(i.e., 


the probability per unit incident flux) for 


detecting a particle C within the phase-space volume element dp\*, A and B are the initial 
colliding particles, C is the produced particle of interest, A' represents all other particles produced 
in the collision, E is the total energy of the produced particle C\ and Q is the solid angle. This 
form (eq. (1)) is favored because the quantity is invariant under Lorentz transformations. 


The data for pion production in proton-proton interactions are primarily reported in terms of 
the kinematic variables, 0 * , y/s, andp^, which are, respectively, the center of momentum (COM) 
frame scattering angle, the invariant mass, and the transverse momentum of the produced pion. 
The mass y/s is a Lorentz invariant quantity and is equal to the total energy in the COM frame; 
p_l_ = p* sin 0* } where p* is the COM momentum; p j_ is invariant under the transformation from 
the laboratory (lab) frame to the COM frame. (See appendix B for a more detailed discussion 
of kinematic variables.) All momenta, energies, and masses are in gigaelectron volts. 


2. 1 . N eut ral P ions 

Biisser et al. (ref. 12) have fitted the LIDCS data obtained in the reaction p + p — tt° + A" 
(where p represents a proton, 7 t g represents the neutral pion produced, and X represents all 
other produced particles) to an equation of the form 


efi <r M _ 
/? = Ap n 

dp° 




exp 


V s 


( 2 ) 


with A = 1.54 x 10“ 2G , n — 8.24, and b — 26,1. Equation (2) is based on a specific set, of 
experimental data, with all measurements taken at 0 * % 90° and was originally intended only 
for pions with high . Comparison of this parameterization with data available from other 
experiments (refs. 13 to 19) indicates that the global behavior of the invariant cross section 
cannot he represented by a function of this form. (See figs. 1 to 96.) The parameterization of 
Biisser et al. (ref. 12) was not plotted because the cross section Is much too small in the pj_ 
range covered by the graph. 

The following form has been used by Albrecht et al. (ref. 20) to represent neutral pion 
production: 

cf_2»_y 

W \PJL + PoJ 
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( 3 ) 


dp* \PJL + PQ, 

where C\ n } and p^ are free parameters. Because this equation only has dependence on p± when 
the data (refs. 13 to 19) shown in figures 1 to 96 also have dependence on \fs and 0 *, this form 
is not general enough to represent all the data. 

Many authors (e.g., ref. 21) have favored a representation for the invariant cross section of 
the form 

17^4 = A (p1 + M*y N/2 f{x ± ,0*) 0 ) 


where f(x j_, 0*) = (1 — x±) J 
Pmax \/s 


N and F are free parameters, the scaling variable x± is given by 




and 


Pmax 


(* + 


to : 


4 iri) 
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where m* and m p are the mass of the neutral pion and the proton, respectively. (See appendix A 
for details on p* nax .) The outline of this basic form has been used by Carey et al. in fitting 
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the invariant cross section for the inclusive reaction p + p — * 7r° + A". (See ref. 22.) Their 
representation is given by 


B |r=''(>>i+ 0 -86) “'“(i— ■*«) 


- 4.5 


dp 


( 5 ) 


where — is the radial scaling variable and the normalization constant A has been 

Pi a ax 

determined as .4 ~ 5. This parameterization accurately reproduces the data for measurements 
taken at 0 * — 90° and yfs > 9.8 GeV but does not agree well with the data for smaller angles 
an d y/s = 7 GeV as can be seen in figures 1 to 96. 

Another problem with this parameterization becomes apparent, when one considers that 
integration over all allowed angles and outgoing particle momenta should yield the total inclusive 
cross section. The details of this calculation appear in section 3. A comparison of the 
experimentally determined total cross section data from Whitmore (ref. 23) with the results of the 
numerical integration of equation (5) shows that the total cross section is greatly underestimated 
by Carey. (See fig. 5.) 

Stephens and Badhwar (ref 19) obtained data from the photon cross sections given by 
Fidecaro et al. (ref. 13), which were taken at incident proton kinetic energy of T\ & ) } = 23 GeV 
and — 0.1— 1.0 GeV. (Note; No error was listed by Fidecaro et al. for pion production. Error 
bars of 10 percent were added to the data in the figures because this level of error was standard 
for most other data. Also, Stephens and Badhwar use the notation E p instead of 7i a p.) Figures 1 
to 96 demonstrate the accuracy of the parameterization of Stephens and Badhwar in this region 
as well as oilier regions. The parameterization of the 7r° invariant cross section proposed by 
Stephens and Badhwar (ref. 19) is presented as follows: 


where 


dp* 


= Af(T My )( l-£)*exp| 


Bp±_ 

1 + Anifjs 



( 6 ) 


Ci-C 2 p± + Qp'1 

\J 1 +4m$/s 



and .4 = 140, B = 5.43, C'i = 6.1 ,C 2 = 3.3, C 3 = 0.6 with xf. = -P-, arK ] = P *cos6*. 

P max II 

The Stephens-Badliwar parameterization was found to be the best of the previously mentioned 
representations because it accurately reproduces the data in the low pj_ region, where the cross 
section is greatest (figs, 1 to 4), and its integration yields accurate values for the total cross 
section (fig. 5). Equation (6) is, however, a poor tool for predicting values of the invariant 
cross section for pj_ > 2 GeV because the value predicted underestimates experimental data by 
~ 1 0 orders of magnitude. (See figs. 6 and 7.) 

No parameterization currently exists that accurately fits the global behavior of the LIDCS 
data. Previous equations have suffered from being too specific to a particular set of experimental 
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data or from failing to reproduce the total cross section upon integration. For these reasons, a 
new parameterization is desired — one that correctly predicts all available data while maintaining 
the essential quality of correctly producing the total cross section upon integration. 

The approach adopted in the present work is to assume the following form for the invariant 
cross section: 

= (sin 0*) D ('/*’P±' 9 *)F(y/7, P± ft*= 90°) (7) 

dp 4 

The motivation for an equation of this form is that as the angle decreases, the cross section 
decreases very slowly at lower values. The approximation that was made in deriving 
equation ( 7 ) is that as pj_ — 0 , the cross section is assumed to be independent of the angle. 

Under the assumption that the invariant cross section can be fitted by equation (7), the 
program goes as follows. Find a representation for the cross section as a function of energy 
\fs and transverse .momentum p±_ from experiment al data taken at 0 * — 90°. The quantity 
F( ^,P±) is then completely determined because (sin#*)^ is unity at 0* — 90°. 

At 0* = 90°, the data are well represented by 

^=90°)= F{^, P± ) ( 8 ) 

dp° 

with 

F(v/s, />_!_) = In (-?=—) G(q.p±) 

V V s ruin / 

7 = * 1/4 

and the COM pion production threshold energy 

in = 2 m v + 777 j,- 


The function 


C(q,p ± ) = 


e4$(0*= 90°) 

dp* 


i„ (-#-) 

V \A s 'nim / 

was parameterized as 

G(q, ]>±) = ex|> h'l + k’2P± + k'i + k 4P l + k$q~' 2 + hp±q~ l + kjp^ + kgq 
+ h p±q~ 2 + k 10 P±q~ l + k 1 lPj 3 ) 


-3 


(9) 


with k i = 3.24, k-2 = —6.046, k$ — 4.35, £4 = 0.883, k~ } = —4.08, kg — —3.05, k 7 = —0.0347, 
kg : 3.046, kg = 4.098, A- j 0 = —1.152, and k\\ = -0.0005. The parameters k\ to were 
obtained with the numerical curve-fitting software, Table Curve 3D (ref. 24), and t lie eleventh 
term was added to modify the low pj_ behavior of the parameterization. 

With F(y/s,p±) determined, the function D (\/s,p±,0*) is the only remaining unknown. 
Solving for D yields 


D{^,P±f*) 


]n^E^J-]n[F(^.p L )} 

1 n (sin 0 ¥ ) 


( 10 ) 
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Equations (8) and (9) were then used in equation (10) to calculate values of D(\fs } j) 

If the function D is independent of angle, then equation (10) could be determined for any fixed 
angle, 0* 90°. Data were compared for a range of angular values, and these data revealed 

that the function D is not independent of angle. The angular dependence turned out to be of 
the form (sin 0*)” 0 * 45 , and 


D(v/W*) = (sin#*) 


— 0.45 


TPj 2 ( \/^) C3 + c 4*^= + “/j + 


1 . 0 ] 


( 11 ) 


with c i = 205.7, c-2 — 3.308, C3 = —2.875, C4 = 10.43, and C5 = 0.8. The final form of our 
resultant parameterization for the neutral pion invariant cross section in proton-proton collisions 
is equation (7) with given in equation (11), F(pj_, v /s) given in equation (8), and 

(7 (r/ given in equation (9). This form is accurate over a much greater range of transverse 
momentum values than those covered by previous representations. (See figs, 1 to 96 for 
comparisons.) For the low transverse momentum region wdiere the cross section is the greatest, 
the fit is quite similar to that of Stephens and Badhwar (ref. 19). Also, figure 5 shows that 
both formulas (eqs. (7) and (6)) integrate to approximately the same total cross section, which 
is in agreement with the data from Whitmore (ref. 23). A more complete comparison of the 
integrated total cross section to data is given by Stephens and Badhwar (ref. 19). 

2.2. Charged Pions 

The available data for charged pions, consisting mostly of measurements made at 0* = 90°, 
are less extensive than 7r° data. Therefore a higher degree of uncertainty exists in LIDCS ’s 
for charged pions. Integration of an LI DCS to get a total cross sec tion and comparison of the 
results with total cross section dat a, allow a check of the global fit of a parameterization. This 
check was made for charged as w r ell as neutral pions, but because of a lack of data, it is more 
important for charged pions. Parameterizations that do not integrate to the correct total cross 
section can be ruled out, even if the LI DCS data are well represented becanse the global behavior 
of the parameterization cannot be accurate. However, producing a correct total cross section 
upon integration does not necessarily imply that the global behavior of the parameterization 
is correct. If more measurements were made, a tighter constraint could be placed on possible 
LIDCS parameterizations. If the spectral distribution is measured at three different values 
of pion energy for two different proton collision energies, the general behavior of the spect ral 
distribution could be checked. The angular dependence of LIDCS parameterization could then 
be tested by integrating over the angle and comparing the results with the spectral distribution 
data. For the purposes of space radiation shielding, measurements at proton lab kinetic energies 
of 3 and 6 GeV and pion lab kinetic energies of 0.01, 0.1, and 1 GeV would be useful because 
this is the region with both a large cross section and large galactic cosmic ray fluxes. These 
measurements would need to be made' only for one pion, preferably 7r°, because the general 
behavior of all the pion production cross sections is approximately the same. With these facts 
in mind, a comparison of LIDCS parameterizations with data from references 15 and 25 to 28 
for charged pion production follows. 

A parameterization for of the form 



( 12 ) 


- P 


has been given by A 1 brow et al. (ref. 28), where A and R are tabulated functions of 

ema 

and A and B are given only for x — 0.18, 0.21, and 0.25, which limits the usefulness of th 
p ar ame t er iz ati on . 
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Alper et. al. (ref. 25) have fitted the data for both tt + and i r production to the following 
form: 

= A exp (-flp . l + C P {) exp {-By 2 ) ( 13) 

where 7 / is the longitudinal rapidity, and A , B ) C } and D are tabulated functions of .s that are 
also dependent on the type of produced particle (7r + or tt~). (Note that at = 90°, y equals 0.) 
The fit to the data is excellent for low transverse momentum, as can be seen in figures 97 to 102, 
but the figures show that this form has an increasing cross section for high pj_, which contradicts 
the trend in the data. Also, there are different sets of constants for each different energy, which 
makes a generalization to arbitrary energies difficult . 


Parameterizations done by Carey et ah (ref. 29) and Ellis and Stroynowski (ref. 21) have a 
similar form, although Carey s was applied only to 7r“. Both underestimate LIDCS's for low p±, 
where the cross section is the largest. (See figs. 97 to 102.) The following equation is Carey s 
p ar ameter iz ati on : 

E ^-j( 7T — ) = N (A + 0.8g) )' 5 ( H) 

where N — 13 is the overall normalization constant and xr = Hp- — & 

Pniax 

equation is Ellis’s parameterization, which was applied to both 7 and 7 r~ production at 
0* = 90°: 

e 7A = yl (d + m2 ) ~ A/2(1 ~ X±)F (15) 


2p 


The following 


where A/, A r , and / are given constants, A is an unspecified overall normalization for. which we 

used A = 13, and ii = ^ ^ ^ jr -. 

Pinax v 5 


The most successful previously developed LIDCS parameterization available for charged pion 
production was found to be the one developed by Badhwar, Stephens, and Golden (ref. 30), 
which is 


F— - zDl r - B >‘ l/(l+4"<p/ s ) 

f/p-l (1 + im'fjs) r 


( 16 ) 


where q is a function of p± and s, sucli that 


Cl + C-2P1 + Cwl 
(1 + Amps) 


and 


*\\ 2+ 7(p± + m ^) 



Here x* = -JS- « 2-31. For tt+ A = 153, B = 5.55, Ci = 5.3667, C 2 = -3.5, C 3 = 0.8334, and 

'■ I 7 max VS 

r = 1. For 7T , A = 127, B = 5.3, Ci = 7.0334, = -4.5, C 3 = 1.667, and r = 3. This form 

is accurate for low transverse momentum (figs. 97 to 102), which is the most important region 
for radiation shielding because of the large cross section. It is also in closed form so that extra 
numerical complexities do not have to be considered. A comparison with a few data, points, 
shown in figures 103 and 104, demonstrates that it integrates to the correct total cross section. 
A more detailed comparison of the integrated cross section with experimental data is given by 
Badhwar, Stephens, and Golden (ref. 30). Because of its relative accuracy and simplicity, this 
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parameterization was integrated to get total cross sections and spectral distributions for charged 
pions. 


Mokhov and Striganov (ref. 31) have also developed the following formulas for both w + and 
7r production: 

where 


Vi = (1 - D) exp(— .EpjJ + Z)exp(-Fpj_) 

v _ 0.2625 

{p\ + 0.87) 


(p± < 0.933 GeV) 


(p ± > 0.933 GeV) 


an d 


V - 1 = 0.7363 exp(0.875pjJ 
Vi = 1 


(P.L < 0.35 GeV) 
(p j > 0.35 GeV) 


with A = 60.1, B = 1.9, and C = 0.18 for ;r + ; ,4 = 51.2, B = 2.6, and C = 0.17 for and 
D = 0.3, E = 12, and F = 2.7 for both tt+ and n~. Figures 97 to 102 show that the formula of 
Badhwar lias a. better fit to the data in the low pj_ region where the cross section is the largest. 

3. Spectral Distributions and Total Cross Sections 

3.1. Method of Generating Other Cross Sections From LIDCS's 

Although L1DCS 's contain all the necessary information for a particular process, sometimes 
other cross sections are needed. For example, one-dimensional radiation transport requires 
probability density distributions that are integrated over solid angle. These quantities are 
calculated in terms of spectral dist ributions and total cross sections rather than LIDCS’s, but 
with accurate parame te rizat io ns of LIDCS’s, formulas for both spectral distributions and total 
cross sections can be developed. LIDCS’s for inclusive pion production in proton-proton collisions 
contain dependence on the energy of the colliding protons y/s, on the energy of the produced 
pion TV, and on the scattering angle of the pion 0. Total cross sections a, which depend onlv 
on y/s, and spectral distributions which depend on s/s and TV, can be extracted from an 

LIDCS by integration. If azimuthal symmetry is assumed, these cross sections take the following 
forms: 


a = 


(1(7 

max 

<7’ ! a 

dE =2 ’ v . 

/ do 
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^ dp^ 

f^n lax 


</V 

7 r / dO 
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/ dp 

^Pniin 

E rfp 3 " 
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(19) 


where 0 max , p ma x> and p [n j n are the extrema of the scattering angle and momentum of t he pion, 
and m w is the rest mass of the pion. 

Tn tl ie COM frame, these extrema can easily be determined. (See appendix A for a detailed 
analysis.) Using conservation of momentum and energy, one can easily show 


> 2 _ {s + rnj-s*) 2 
45 


m 7 


( 20 ) 
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where s x is the square of the invariant mass of the sum of all particles excluding the pion, and p 
is the magnitude of the three momentum of the pion. The independence of p on 0 implies that 
6 can take on all possible values (i.e $ IIiax = 7r), and the symmetry of the COM frame implies 
that p l[lin = 0. For a given value of s y it is obvious that moment um is a maximum when s x is 
a minimum. As shown in appendix B, an invariant mass is a minimum when it is equal to the 
square of the sum of the rest masses of the particles in question. Momentum is, therefore, a 
maximum when s x is the square of the sum of the least massive combination of particles that can 
be produced and still satisfy all relevant, conservation laws. For the reaction p + p — *■ 7r + .r, we 
h av e .s' j- ^ 47772 , where subscript p represents a, proton. Exact formulas are listed in appendix B 
along with a more detailed analysis. 

If a Lorentz transformation is applied to the maximum COM momentum, the integration 
limits can be determined in other frames. Byckling and Kajantie (ref. 32) have shown that by 
transforming to the lab frame, the following formula can be obtained: 




Pa ®max 


Vs cos 0 ± (E a + m p )\f 


*Pmax 


m nPa Si 11^ 0 


$+ pl sin 2 (0) 


-1 


( 21 ) 


where quantities with an asterisk are COM variables, quantities without an asterisk are either 
lab or invariant variables, rv^ Is the rest mass of a proton, p (l is the magnitude of the moment um 
of the projectile proton, and p+ ~ p max Is the maximum pion moment um. The greater of the 
two quantities 7 ?” = p m j Tl and 0 is the minimum pion momentum, and the maximum scattering 
angle can be determined by the requirement that p± be real. This requirement implies that the 
quantity under the square root in equation (21) must be greater than or equal to 0. Solving for 
0 , nax then gives the formula 

fl.nax = sill " 1 ( ^ ) ( 22 ) 

V Pa m V ) 

With the limits of integration determined, an LIDCS can be turned into a total cross section 
or a spectral distribution by numerical integration. This procedure will, however, give discrete 
“data” points not closed-form expressions. Parameterizations of these numerical data are needed, 
if relatively simple formulas for these cross sections are desired. This process was completed for 
all three pion species, and the corresponding formulas are given in the next section. 

3.2. Parameterizations 

The surface parameterizations for the spectral distribution as a function of incident proton 
kinetic energy in the lab frame 7]^, and the lab kinetic energy of the produced pion T n have 
been completed by numerically integrating LIDCS charged pion parameterizations of Badhwar, 
Stephens, and Golden (eq. (1 6) and ref. 30) and the neutral pion cross section from Stephens and 
Badhwar (eq. ( 6 ) and ref. 19). The numerical integration routines were checked by computing 
total cross sections in both the lab and COM frames and comparing the results. Because total 
cross section is a Lorentz invariant, the results should be the same in both frames. To accurately 
fit the integration points for low energies, considering two regions of the surface and determining 
representations for them individually have been necessaiy. For each of t he three pions, the two 
regions consist of laboratory kinetic energies 7j a j, from 0.3 to 2 GeV and from 2 to 50 GeV. 
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The neutral pion spectral distribution for the range 0.3 to 2 Gev is given by the following 
equations: 

F 2 = A l T?* + A 3 T** 

F l = ex I> + ^9 lu + A 1 1 Ti A 1 > ( 23) 

= Mlj^+Au exp(.4i CV ^V + .4i7^ 4,8 T 1 ^ 9 ) T^ i:> 

with constants A; given in table 1. The neutral pion spectral distribution for the range 2 to 
50 GeV is given by the following equations: 


F 2 = B<i¥* + B,Tfj 

F[ = ex P ( 4 7 = 


7== + B ~ T & + b 9 t - w + B n^ n 




\dEJ , ab O 

with constants /?; given in table 2. 

Table 1. Const ants ,4, for Equations (23) 


Table 2. Constants /?, for Equations (24) 


Constant 

Valu e 

A\ 

6.78 x IO -10 

M 

-2.86 

A-i 

1.82 x 1(T 8 

A , 

-1.92 

Ai 

22.3 

A, 

0.226 

-4 7 

-0.33 

-4s 

-1.75 

-Ay 

-32.1 

-4 io 

0.0938 


-23.7 

-4 12 

0.0313 

.4 13 

2.5 x 10 fi 

-4 1 1 

1.38 

-4 15 

0.25 

A 16 

-39.4 

-4 17 

2.88 

•4 is 

0.025 

4|y 

0.75 


Constant 

Value 

B i 

1.3 x 10“ 10 


-2.86 

B 3 

4.27 x 10 -9 

B-\ 

-2.4 

B:> 

22.3 

B 6 

-1.87 

B 7 

1.28 

B s 

- 1 .25 

By 

-33.2 

flio 

0.0938 

flu 

-23.6 

fll2 

0.0313 

A 13 

2.5 x 10 6 

Am 

0.25 

fli j 

60322 

fll6 

1.07 

A 17 

-67.5 


I 











The positively charged pi on spectral distribution for the range 0-3 to 2 GeV is given by the 
following equations: 


f-2 = cyr^ + c^ 


F i = exp 


C 5 + -j= + c 7 t,2 + C 9 T?' 0 + Ci 1 7f + Ci 4 In (r lab ) 
v^lab 


> (25) 


(£) =C a lf<°!± + C, 7 7fex, 

V /L /kb *2 


with constants C\ given in table 3. 

The positively charged pion spectral distribution for the range 2 to 50 GeV is given by the 
following equations: 


F-2 = DiT?* + D»7] 


F\ = exp Drj + 


Da 


G + DjT?* + M)7f 10 


= Di\tf' 2 J- + D ls T?'< exp(z? 15V ^V + 0 16 T£j“) 


lal) ^2 

with constants /?; given in table 4. 


(26) 


Table 3. Constants Cj for Equations (25) 


Constant 

Valu e 

C, 

2.2 x 1(T 8 

C-2 

-2.7 

C 3 

4.22 x 10“ 7 

C.1 

-1.88 

C 5 

22.3 

Ce 

1.98 

Ci 

-0.28 

c 8 

-1.75 

c 9 

-29.4 

C, 0 

0.0938 

c„ 

-24.4 

C,2 

0.0312 

Cl3 

0.0389 

C,,, 

1.78 

Cl 5 

2.5 x 10 6 

C, 6 

0.25 

C, 7 

976 

Cl 8 

2.3 

Cl 9 

-46 

C 20 

-0.989 


Table 4. Constants Dj for Equations (26) 


Constant 

Value 

D\ 

4.5 x 10 _u 

D 2 

-2.98 

Ih 

1.18 x 10~ 9 

£> , 

-2.55 

Do 

22.3 

Da 

JO 

o 

1 

D 7 

-35.3 

Ds 

0.0938 

A) 

-22.5 

D iu 

0.0313 

£>n 

2.5 x 10 6 

£>12 

0.25 

A 3 

60322 

£>i i 

1.18 

£>15 

-72.2 

£>16 

0.941 

£>17 

0.1 
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The negatively charged pion spectral distribution for the range 0.3 to 2 GeV is given by the 
following equations: 


F 2 = C 1 T^+G 3 T^ j 

F\ = exp ^Grj + — ■■ + G-T^’ s -f G'c|T£' iu ^ 


(27) 


dEJ lab 


) =7f» 


G]2~ + G43 exp^6’i4 s/lx'j 


with constants (7/ given in table 5. 

The negatively charged pion spectral distribution for the range 2 to 50 GeV Is given by the 
following equations: 


F 2 = FI 

F, = expf H Ti + JjjL. + H 7 T?« + // 9 t " 10 j 

with constants /7; given in table 6. 


Table 5. Constants Gj for Equations (27) 


Constant 

Value 

G i 

1.06 x 10~ 9 

g 2 

- 2.8 

Q 

3.7 x 10 -8 

G\ 

-1.89 

G 5 

22.3 

Q > 

-1.5 

G 7 

-30.5 

g 8 

0.0938 

G 9 

-24.6 

Gio 

0.0313 

Gn 

0.25 

G] 2 

2.5 x 10° 

G 13 

7.96 

Gi i 

-49.5 


Table 6 . Const ants II j for Equations (28) 


Constant 

V alue 

n i 

2.39 x 10 ~ 10 

n-i 

-2.8 

Hi 

1.14 x 10“ 8 

Hi 

-2.3 

Hr, 

22.3 

Ik 

-2.23 

Ik 

-31.3 

Ik 

0.0938 

H g 

-24.9 

H io 

0.0313 

II n 

2.5 x 10 6 

Hn 

0.025 

II 13 

60322 

//n 

1.1 

7/ 15 

-65.9 

77,6 

-9.39 

77 1 7 

-1.25 
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Tbtal inclusive cross sections are given by the following equations: 


/x^q — 0.00/ “l - 0.1 


ln(T lal) ) , 0£ 

Tl 


-l 


+ 


' lab 


IhI> 


( 29 ) 


= 0.00717 + 0.0652 


MZU) , 0.162 


-1 


T lal> 


+ 


T 2 

1 lal> 


(30) 


o- T - = 0.00456 + 


0.0846 0.577 


-l 


1 lab 


+ 


j 1 1 . 5 
Mai) 


(31) 


For neutral pions, spectral distributions and total cross sections that, were based on tin' 
present, parameterizat ion gi von in equation (7) were also developed. The formula for the spectral 
distribution was not divided into two regions and is much simpler than the previous formulas: 


da 

dE 


= exp 

lal> 



K-2 






(32) 


whore I\\ = —5.8, K 2 = — 1.82, A3 = 13.5, and A.j = —4.5. 

Because equation (7) and Stephens' LIDCS parameterization integrate to nearly the same 
total cross section (fig. 5), separate total cross section parameters ations are not necessary (i.e., 
use eq . (29)). 


3.3. Discussion of Spectral Distributions and Total Cross Sections 


As discussed previously, figures 1 to 4 and 6 to 9 show the LIDCS parameterizations for tt° 
production of Carey et al. (eq. (5) and ref. 14), Stephens and Badhwar (eq. (6) and ref. 19), 
and of equation (7) plotted with data from references 12 to 17, 20, and 23. The figures are 


graphs of cross .section 



plotted against transverse momentum p± for various values of 


COM energy E cm and COM scattering angle 0*, which can be transformed into lab variables 
as shown in appendix B. Figures 1 to 3 and 7 show that the parameterization of Carey et al. 
is not a.11 adequate representation of the data. Figures 6 and 7 show that t he parameterization 
of Stephens and Badhwar fails for high transverse momentum by severely under predicting the 
cross section. 


Figure 5 shows numerically integrated LIDCS parameterizations of Stephens and Badhwar 
(eq. (6) and ref. 19), of Carey et al. (eq. (5) and ref. 14), and of equation (7) (referred to 
as “Kruger”) for 7r° production plotted with a parameterization of the integrated formulas of 
Stephens and Badhwar, referred to as “Stephens-total param” (eq. (29)). Three data points 
from Whitmore (ref. 23) show that Carey s parameterization does not integrate to the correct 
values and that the rest are quite accurate. (See ref. 19 for more detail.) 

Figures 10 to 96 show tt 0 LIDCS parameterizations of Carey et al. (eq. (5) and ref. 14), 
of Stephens and Badhwar (eq. (6) and ref. 19), and of equation (7) plotted with data from 
references 12 to 17, 19, and 20 over a wide range of angles and energies. These graphs show that 
equation (7) has the best global fit to all data available. 

As discussed previously, figures 97 to 102 show tt+ and LIDCS parameterizat ions of Alper 
et al. (eq. (13) and ref. 25), of Badhwar, Stephens, and Golden (eq. (16) and ref. 30), of Ellis 
and Stroynowski (eq. (15) and ref. 21), of Carey et al. (eq. (14) and ref. 29), and of Mokliov 
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and Striganov (eq. (17) and ref. 31) and LIDCS dat a from references 18 and 15 plotted against 
transverse momentum for different values of COM energy E cm with all at 0* = 90°. These 
graphs show that the parameterizations of Badhwar and Alper best fit the data, but Alpers 
parameterization rapidly increases for high transverse momentum which contradicts the trend 
of the data. 


Figure 103 shows numerically integrated LIDCS parameterizations of Badhwar, Stephens, and 
Golden (eq. (16) and ref, 30) for tt + plotted with parameterizations of the integrated formulas 
of Badhwar referred to as “present work" (eq. (30)). Figure 104 shows numerically integrated 
LIDCS parameterizations of Badhwar, Stephens, and Golden (eq. (16) and ref. 30), of Carey 
et al. (eq. (14) and ref. 29) for ir~ plotted with parameterizations of the integrated formulas of 
Badhwar referred to as “present work" (eq. (31)). Three data points from reference 23 show that 
Carey’s parameterization does not integrate to the correct values and that Badh war’s formula is 
fairly accurate. The figures also show that the parameterization fits the numerically integrated 
formulas very well. 


Figure 105 shows tt^ spectral distribution parameterizations given by equations (23) and (24) 
plotted with LIDCS parameterization of Stephens numerically integrated at several lab kinetic 
energies. Figure 106 is the same as figure 105 except that the spectral distribution of 
equation (32) is plotted with the numerical integration of equation (7). 


Figures 107 and 108 show tt and 7r + spectral distribution parameterizations plotted with 
LIDCS parameterization of Badhwar, Stephens, and Golden (eq. (16) and ref. 30) numerically 

integrated. Cross section is plotted against the kinetic energy of the produced pion T tt at 

several values of the lab kinetic energies of the colliding proton. 


4. Concluding Remarks 

This paper presents parameterization of cross sections for inclusive pion production in proton- 
proton collisions. The cross sections of interest are Lorentz invariant differential cross sections 
(LIDCS's), laboratory (lab) frame spectral distributions, and total cross sections. For neutral 
pious the parameterization of Stephens and Badhwar (Asiwphys. & .Spare Sci. y vol. 76, 1981, 
pp. 213 217) fit the data well for low values of transverse momentum pj_ but overpredict the 
cross section by many orders of magnitude at high p± values. Because of this inaccuracy, an 
equation was developed. The final form of our resultant parameterization for the neutral pion 
invariant cross section in proton-proton collisions is as accurate as that of Stephens and Badhwar 
at low p_L values but is much more accurate at high pj_ values. For charged pions the formula of 
Badhwar, Stephens, and Golden (Phys. Rev . D ) vol. 15, 1977, pp. 820 831) was found to best 
represent, the data except at high p± values and that of Ellis and Stroynowski (Rev. Modem 
Phys vol. 49, 1977, pp. 753 775) was quite accurate. The formula of Badhwar, Stephens, and 
Golden was used in the development of spectral distributions and total cross sections because it 
was the most accurate at low pj_ where the cross section is the greatest. 

The data for lab frame spect ral distributions and total cross sections are scarce; therefore, 
parameterizations for these quantities were developed with LIDCS formulas. Those formulas 
were numerically integrated, resulting in discrete numerical data points for the other cross sec- 
tions, namely spectral distributions and total cross sections. The accuracy of the representations 
of lab frame spectral distributions and total cross sections is, therefore, limited to the accu- 
racy of the original LIDCS’s. The numerical data were then parameterized so that closed-form 
expressions could be obtained. As a check on the accuracy, the total cross section numerical 
data were compared with experimental data. They were found to agree quite well, but when 
the numerical data for the spectral distributions for the formulas for 7r° production are com- 
pared, they are found to disagree. Because both original LIDCS formulas fit the data well at 
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low p _ |_, where the cross section is greatest, and both formulas integrate to the correct total cross 
section, the available data must not be sufficient to uniquely determine the global behavior of 
the LIDCS’s. The data for charged pion production weie much more limited than the data for 
neutral pion production; therefore, the same problem exists for charged pions. 

To more accurately determine the cross sections for space radiation applications, measure- 
ments of the spectral distribution at lower energies (for example, proton lab kinetic energies of 3 
and 8 GeV and pion lab kinetic energies of 0.01, 0.1, and 1 GeV) for one pion species would 
need to be taken. Only measurements for one pion species would be needed because they all 
have approximately the same general behavior. These measurements would put a much t ighter 
constraint on the global properties of the LIDOS 's, and the spectral distribution parameteriza- 
tion s could also be made more accurate. 
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Appendix A 

Derivation of Maximum Momentum 

Consider an inelastic two-particle collision. The maximum momentum that a produced 
particle can have in the COM frame can easily be determined by imposing conservation laws. 
The reaction considered is a two-particle A and B reaction resulting in a pion tt and various 
other particles X (i.e., A + B — >■ n + X). Quantities pertaining to the initial particles are labeled 
with subscripts a and b. Quantities pertaining to the pion of interest are labeled with subscript 
7T, and quantities pertaining to the system consisting of all other produced particles are labeled 
with subscript x. Units where the speed of light is equal to unity are used. 

Conservation of energy implies that the initial energy equals the final energy: 

E a + Eb = E z + E x ~ y/s ( A 1 ) 

where \fs is the invariant mass of the entire system. Sec appendix B for the relations between 
various kinematic variables. 

Tn the COM frame, the total three moment um is zero. Therefore, 


J3i 

+ 

II 

Pit + Vj = 0 

(A2) 

Pit = 

-Pr 

(A3) 

|Px| 2 = 

\pA 2 

(A4) 

El ~ ml = 

El - s r 

(A5) 


where s x is the square of the invariant mass of the system consisting of all particles except the 


pion, and 777 ^ is the rest mass of the pion. By rearranging terms, adding 
a little further algebra, equation (37) becomes 

E'l to both sides, and 

E*+E* = 2E* + a r -ml 

(A6) 

El + El + 2 E* E x =2 E* + k x - ml + 2 E X E X 

(AT) 

s Sj- + THjj. — 2E-x{E n -f- E x ) 

(AS) 

Substituting equation (A8) into equation (Al) results in 


a«ar O 

2 E^s/s = s — s x - h m £ 

(A9) 

9 

^ S + 777 t. — S T 


E 

(A10) 

P ,| 2 = <S + " S ' )2 ml 

(All) 


Equation (All) obviously implies that the pion momentum is a maximum for a given s, w r hen 
s x is a minimum: 


IP*n 


4s 




(A12) 


16 


To specify |p 7 r n i ax p for a given inclusive reaction, needs to be derived; s r is the square of 
the invariant mass of an N particle system, which is defined as the square of the four momentum 
of that system, as follows: 



sx = 2 ! Pi ■ pj 
0=1 


(A 13) 


(A 14) 


where p; pj is a Lorentz scalar and can be calculated in any frame. For simplicity, the calculation 
is done in the rest frame of the ith particle, so that Ej = mp the speed v / = 0, p; = 0, and 


V - 



For massive particles, 


where 7 j 


Pi ■ Pj = E i E j - Pi ■ Pj (A15) 

Pi pj = imnjmj (A 16) 

Pi ■ Pj = lj m i m j (A 17) 



> 1 , and vj is the speed of the jth particle. Equation (A17) further implies 


Pi * Pj > mpm j 


(A18) 


If one particle is massless, pi pj can be calculated in the rest frame of the massive particle, 
which can be assumed to be the ith particle, without a loss of generality: 

Pi ■ Pj = E i E j - Pi ■ Pj ( A 1 9) 

Pi ■ Pj = rrt iEj ( i / j) (A20) 

which is a minimum and equals zero in the limit as Ej — • 0. If both particles are massless, 


Pi Pj = EiEj ~ Pi ■ Pj ( A‘2 1 ) 

Pi ■ Pj = E i E j ~ Ej Ej ( A 22 ) 

Pi ' Pj — 0 (A23) 

Therefore, for any combination of massive and massless particles, 

Pi ■ Pj > m;mj (A 24) 


Substituting equation (A24) into equation (A 14) results in 



Therefore, 


8 


min 


N 

J2 m < m j 


0=1 


s 


min 



( A 26 ) 


(A 27) 


when the sum is over the least massive combination of particles that can be produced while 
satisfying all relevant conservation laws. 


For pions produced in proton-proton collisions, the reactions where s x ■ can be a minimum are 
as follows: 


For ? r + , 


For 7r^, 


For 7r , 


P + P 7T + + p + p + € + V t 

(A28) 

=> •‘ ! J ln in = i' 2vi v + ™C + ™—) 2 ~ 

(A29) 

p + p— +p + p 

(A30) 

=> Sf . — 4 mjj 

a mi n r 

(A31) 

p + p — ► 7r~ + p+ p-\- e + + 

(A32) 

=» *x min = (2 m p + m t + m lJt ) 2 ss 4?;^ 

(A33) 


where p represents a proton, e represents an electron, represents an electron-neutrino, V t 
represents an elect- ron-antineutrino, and represents a positron. 


18 



Appendix B 
Kinematic Relations 

The two reference frames of interest in this paper are the COM frame and the lab frame. In 
a two-particle collision, the lab frame is defined so that one of the colliding particles is at rest, 
and in the COM frame, the net three momentum is zero. Using the fact that the square of the 
four momentum of a system s is Lorentz invariant , a relation between COM energies and lab 
energies can be derived. In the COM frame, 

s = (E* a + E* b )‘ 2 -(fi + Plf 

= (£,: + Elf ' (Rl) 

where * refers to a COM quantity, and variables without * are either in the lab frame or they are 
invariant under the transformation from the lab frame to the COM frame, E is a total energy, 
p is a three moment um, and subscripts a and b distinguish the two colliding particles. Equation 
(Bl), therefore, implies that \/s is equal to the COM total energy. 

In the lab frame where particle b is at rest, and T| a |, is the kinetic energy of particle a, 


s — (E a + Ejj f — (p tt + p;,) 2 

(B2) 

s = (E a + m b ) 2 - |p„| 2 

(B3) 

s = 7U 2 + 771 2 + 2E a m b 

(Bl) 

S - 7?7 2 + 77 ? 2 + 2(T| a |, + m„)mi 

(B5) 

s - (m a + m b ) 2 + 2T ]ilh m b 

m 

where m„ and mj, are the rest masses of particles a and b, respectively, Tj a j 
the lab frame kinetic energy of the incoming particle, and E = \j\vf + , " 2 - 
equations (Bl) and (B6) imply that 

, = m a is 

Taken together, 


(E* + E* b f = (m„ + m b f + 2T lal) m 6 (B7) 

The relations between the lab scattering angle 0, the magnitude of the lab three momentum 
|p|, and the corresponding COM quantities can he derived by using the following equations 
(ref. 32): 

|p|*sin 0* = |p| sin 0 (B8) 

|p|* cos 6* = 7 |p | cos 0 — 7 v E (B9) 

E* = —7 v |p| cos 0 + 7 E (B 1 0) 
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where v = — 7 2 is the relative speed of the lab frame and the COM frame, and for the case 

T 1 2 jjj 

of two colliding protons, 7 — — — a- — ( ref. 32). The magnitude of the COM three momentum 

V 5 

and the COM scattering angle can now be easily derived in the following manner: 


and 


The refore, 


|p|* 2 = |p|* 2 (cos 2 0* + si»V) 

|p| 2 = 7 2 ( \p\ cos 0 — vEf + |p| 2 sin 2 0 


tan 6* — 
tan 0 — 


|p|*sin 6 * 

|p[* cos 6 * 

|p | sin 0 

7 \p\ cos 6 — jvE 


0 * = tan 


-1 


|p| sin 0 


r |p| cos 0 — 7 vE 


(Bll) 
(B 12) 

(B13) 

(BH) 

(B 15) 
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Appendix C 

Synopsis of Data Transformations 

The data used in the comparison of different parameterizations were given in terms of several 
different kinematic variables. Some data were transformed so that all data would be expressed 
in terms of the same variables. A synopsis of the transformations that were performed for the 
data, plotted in the figures is as follows. 

The data from Biisser et al. (refs. 12 and 15) and Owen et al. (ref. 17) were given for 
different values of >/s, pj_, and 8* . No transformations were performed. The data from Carey 
et al. (ref. 14) were given for different values of Pp, pj_, and 8\ Pp was transformed to \fs 

by using equation (B4) and the relation E — yj |pp + m 2 . Then 8 was transformed to 8 * with 

equation (B 15). Eggert et al. (ref. 10) used y/s, pj_, and 8 . Equation (B 15) was used to transform 
8 to 8 *. Alper et al. (ref. 25) used y/s, pj_, and the longitudinal rapidity y . Only data with 
y — 0 were used, and when y = 0 then 8 * equals 90°. The data from Whitmore (ref. 23) were 
not transformed. 

Stephens and Badhwar (ref. 19) used photon production data from Fidecaro et al. (ref. 13) 
to derive pion production cross sections. The variables T \ d \) . 8 , and p were used by Stephens. 
Equation (B6) was used to transform T\. d \ y into >/s. Next, 0 was transformed to 8 * by using 
equation (B15). Finally, p was transformed into pi by using the equation pi — psmO. 
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Figure 1. LIDCS plotted against transverse momentum for i r° production for E rm = 7 GeV and 6 
with data at 12.2° < 9* < 12.4°. 





Figure 2. LIDCS plotted against transverse momentum for tt° production for E cm — 7 GeV and 0 
with data at 40 .3° < 0* < 4 1.9°. 






Figure 3. LIDCS plotted against transverse momentum for 7r° production for E rm = 7 GeV and 0 
with data at 88.3° < 6* < 90°. 



Figure 4. LIDCS plotted against transverse momentum for 7 r° production for E nv — 9.8 GeV and 0 
with data at 37.1° < 0* < 37.5°. 





Figure 5. Parameterization of total tt° production cross section plotted with numerically integrated LIDCS 


p ar a 1 1 le t e ri z at ion s . 



Figure 6. LIDCS plotted against transverse momentum for tt° production for E crn = 53 GeV and 0* = 90° 
with data at 0 * — 90°. 
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Figure 7. LIDC'S plotted against transverse momentum for tt° production for E cm — 53 GeV and 0 
with data at 0* — 5°. 



Figure 8. LIDCS plotted against transverse momentum for tt 0 production for E nv = 9.8 GeV and 0* 
with data at 86.8° < 0* < 87.9°. 







Figure 10. LIDCS plotted against transverse momentum for tt 0 production for E r m = 12 GeV and 9 
with data at 21.4° < 9* < 21.6°. 








Figure 12.LIDCS plotted against transverse momentum for 5r° production for E nv = 12GeV and $ 
with data at 44.7° < 9* < 45.3°. 
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Figure 15. LTDCS plotted against transverse momentum for 7r° production for E rw — 12 GeV and 0* = 65° 
with data at 64.6° < 9* < 66.2°. 







Figure 18. LIDC'S plotted against transverse momentum for 7r° production for E cm = 14 GeV and#’ 
with data at 40.1° < 0* < 40.6°. 







p v GeV 

Figure 19. LIDCS plotted against transverse momentum for tt° production for E cw = 14 GeV and 9 
with data at 50.7° < 9* < 52.6°. 
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Figure 20. LIDCS plotted against, transverse momentum for 7r° production for E cm — 14 GeV and# 
with data at 60.6° < #* < 62.1°. 







Figure 22. LIDCS plotted against transverse momentum for 7r° production for E cm — 14 GeV and 9 
with data at 72.3° < 9* <74.3°. 







p ± , GeV 

Figure 24. LTDCS plotted against transverse momentum for 7T° production for E cm = 15.4 GeV and 
0* = 28° with data at 27.5° < 0* < 27.7°. 
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Figure 25. LIDCS plotted against transverse momentum for 7r° production for E r w 
0 * = 45° with dat a at 44.4° < d* < 47.7°. 
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Figure 26. LIDCS plotted against transverse momentum for 7r° production for E cm = 15.4 GeV and 
0* rr 56° with data at 55.9° < 0* < 57°. 
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* Experimental data (ref. 14) 

Carey (eq. (5) and ref. 14) 

Stephens (eq. (6) and ref, 19) 

Present work (eq. (7)) 
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Figure 27. LIDCS plotted against transverse momentum for w 11 production for E nll = 15.4 GeV and 
6* = 61° with data at 59.6° < 0* < 62.8°. 



Figure 28. LIDCS plotted against transverse momentum for tt° production for Ecm 
0* = 67° with data at 66.5° < 0* < 68.3°. 
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Figure 29. IJDCS plotted against transverse momentum for 7r n production for E rw = 15.4 GeV and 
0 * — 79° with data at 78.5° < 0* < 80.8°. 



Figure 30. LIDCS plotted against transverse momentum for 7r° production for E C m 
0* — 30° with data at 30.0° < 0* < 30.3°. 
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* Experimental data (ref. 14) 
— Carey (eq. (5) and ref. 14) 

■ — Stephens (eq. (6) and ref. 19) 
Present work (eq. (7)) 
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Figure 31. LTDCS plotted against transverse momentum for 7r° production for E cm = 16.9 GeV and 
r = 48° with data at 42.2° < 6* < 48.9°. 


★ Experimental data (ref. 14) 

— Carey (eq. (5) and ref. 14) 

— Stephens (eq. (6) and ref. 19) 
Present work (eq. (7)) 
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Figure 32. LTDCS plotted against transverse momentum for tt° production for E nv = 16.9 GeV and 
0 * = 61° with data at 55.1° < 0 * <61.6°. 





Figure 33. LIDCS plotted against transverse momentum for 7r° production for E cw = 16.9 GeV and 
0* = 72° with data at 71.3° < 0* < 73.4°. 



Figure 34. LIDCS plotted against transverse momentum for 7r° production for E cm = 16.9 GeV and 
0* = 84° with data at 83.6° < 0* < 86.2°. 
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* Experimental data (ref. 14) 

Carey (eq. (.*>)and ref. 14) 

Stephens (eq. (6) and ref. 19) 

Present work (eq. (7)) 
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Figure 35. LIDCS plotted against transverse momentum for tt^ production for E rw — 17.3 GcV and 
0 * = 84° with data at 81.2° < 0* < 84.3°. 



Figure 36. LIDCS plotted against transverse momentum for 7r° production for E cw 
0* = 32° with data at 32.3° < 0* < 32.5°. 
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Figure 37. LIDCS plotted against transverse momentum for production for Ecw — 18.2 GeV and 
Q* — 40° with data at 39.2° < 0 * <41°. 



Figure 38. LIDCS plotted against transverse momentum for i r° production for E cm = 18.2 GeV and 
6 * = 52° with data at 51 .5° < 0* < 54.6°. 
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* Experimental data (ref. 14) 

— Carey (eq. (5) and ref. 14) 

Stephens (eq. (6) and ref. 19) 

Present work (eq. (7)) 
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Figure 39. LIDCS plotted against transverse momentum for tt° production for E cm — 18.2 GeV and 
0* = 64° with data at 64.2° < 0* < 65.7°. 



Figure 40. LIDCS plotted against transverse momentum for tt 0 production for E rw 
0* = 76° with data at 75.5° < 0* < 77.9°. 
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Figure 41. LIDCS plotted against transverse momentum for x° production for Ecm = 18.2 GeV and 
0* = 88° with data at 88° < 0* < 89.2°. 



Figure 42. LIDCS plotted against transverse momentum for t (i production for Enn = 18.6 GeV and 
0* = 80° with data at 76.9° < 0* < 80.1°. 






Figure 44. LIDCS plotted against transverse momentum for production for E C m = 19.4 GeV and 
6* = 50° with data at 48.4° < 9* < 51.5°. 





Figure 45. LIDCS plotted against transverse momentum for 7r° production for E cw = 19.4 GeV and 
0 * — 55° with data at 54.6° < 0* < 55.6°. 



Figure 46. FIDOS plotted against transverse momentum for 7r° production for E nn 
0* = 68° with data at 67.7° < 8* < 69.4°. 
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Figure 47. LIDCS plotted against transverse momentum for 7r° production for Erm — 19.4 GeY and 
0* — 80° with data at 79.1° < 0* <81.6°. 



Figure 48. LIDCS plotted against transverse momentum for tt 0 production for E rm 
9 * = 87° with data at 85.1° < 0* < 88.1°. 
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Figure 49. LIDCS plotted against transverse momentum for tt° production for Ecm — 19.8 GeV and 
0* = 75° with data at 73.2° < 0* < 76.5°. 



Figure 50. LIDCS plotted against transverse momentum for 7r° production for E nn 
0* — 36° with data at 36.3° < 0* < 36.8°. 
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Figure 51. LIDCS plotted against transverse momentum for 7r° production for Emi = 21.0 GeV and 
0* = 73° witli data at 70° < 0* < 73.4°. 



Figure 52. LIDCS plot ted against transverse momentum for 7r° production for E rm 
0* = 34° with data at 33.2° < 0* < 36.4° and E cm - 20.6 and 21.7 GeV. 


= 21.7 GeV and 
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Figme 53. LTDCS plotted against transverse momentum for tt° production for E cw — 21.7 GeV and 
0* — 45° with data at 43.7° < 0* < 48.9° and E cm — 20.6 and 21.7 GeV. 
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Figure 54. LIDCS plotted against transverse momentum for i r° production for E C m =21.7 GeV and 
0* = 60° with data at 57.4° < 0* < 61.2° and E r ™ = 20.6 and 21.7 GeV. 
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Figure 55. LIDCS plotted against transverse momentum for tt° production for E rw — 21.7 GeV and 
0* = 74° with data at 70.9° < 9* < 75.8° and E cm = 20.6 and 21.7 GeV. 



Figure 56. LIDCS plotted against transverse momentum for r° production for E C m = 21.7 GeV and 
0* = 80° with data at 78.5° < 0* < 85.2° and E m = 20.6 and 21.7 GeV. 





Figure 57. LTDCS plotted against transverse momentum for tt° production for Ecm 
9* = 86° with data at 85.5° < 0* < 88.4°. 
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Figure 58. LIDCS plotted against transverse momentum for tt 0 production for E rm 
0 * rr 70° with data at 67.2° < 0* < 70.7°. 
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Figure 01. LTDCS plotted against transverse momentum for 7r° production for E cm — 22.8 GeV and 
0 * = 44° with data at 41 .8° < 0 * < 44.8°. 



Figure 62. LIDCS plotted against transverse momentum for tt° production for Ecm =22.8 GeV and 
0 * = 63° with data at 62.4° < 0 * < 63.7°. 
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* Experimental data (ref. 14) ^ 

Carey (eq. (5) and ref. 14) 

Stephens (eq. (6) and ref. 19) 

Present work (eq. (7)) 
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Figure 63. LIDCS plotted against transverse momentum for 7r° production for E nv = 22.8 GeV 
0 * = 77° with data at 76.4° < 0* < 78.5°. 



Figure 64. LIDCS plotted against transverse momentum for 7 r° production for Enn ~ 22.8 GeV and 
0* = 89° with data at 88.2° < 9* < 88.9°. 
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Figure 65. LTDCS plotted against transverse momentum for 7r° production for E nn = 23.1 GeV and 
0 * = 67° w itli data at 64.7° < 0* < 68.2°. 



Figure 66. LIDCS plotted against transverse momentum for ? r° production for E cm 
0* — 16° with data at 6* = 15° and 17.5°. 
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Figure 67. LIDCS plotted against transverse momentum for 7T° production for Erw — 23.3 GeV and 
0* = 21° with data at 0 * = 20° and 22°. 



Figure 68. LIDCS plotted against transverse momentum for i r° production for E cm 
0* = 90° with data at 86.1° < 0* < 90° and 23.5 < Em, < 23.8 GeV. 
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* Experimental data (ref. 14) 

— Carey (eq. (5) and ref. 14) 

— Stephens (eq. (6) and ref. 19) 

— Present work (eq. (7)) 
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Figure 69. LIDCS plotted against transverse momentum for 7r° production for E cm — 23.8 GeV and 
0* = 31° with data at 30.4° < 0* < 31.9°. 



Figure 70. LIDCS plotted against transverse momentum for i r° production for E rm = 23.8 GeV and 
6 * = 42° with data at 41.5° < 0* < 42.1°. 
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Figure 71. LIDCS plotted against transverse momentum for tt 0 production for E cm — 23.8 GeV and 
0* — 65° with data at 54 .6° < 0* < 66.0°. 
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Figure 72. LIDCS plotted against transverse momentum for 7r° production for E rm — 23.8 GeV and 
e* = 66° with data at 63.2° < 0* < 66.7°. 
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* Experimental data (ref. 14) 

Carey (eq. (5) and ref. 14) 

Stephens (eq. (6) and ref. 19) 

Present work (eq. (7)) 
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Figure 73. LIDCS plotted against transverse momentum for tt 0 production for E rm — 23.8 GeV and 
9* — 76° with dat a at 73.2° < 0* < 76.6°. 



Figure 74. LIDCS plotted against transverse momentum for 7 r° production for E r ,„ 
9* = 73° with data at 71 .0° < 9* < 74.4°. 
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* Experimental data (ref. 14) 
Carey (eq. (5) and ref. 14) 


Stephens (eq. (6) and ref. 19) 

Present work (eq. (7)) 
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Figure 75. LTDCS plotted against transverse momentum for 7r° production for E nv — 25.7 GeV and 
0* = 71° with data at 68.8° < 0* < 72.4°. 



Figure 76. LIDCS plotted against transverse momentum for 7r° production for E nv = 26.6 GeV and 
0 * = 70° with data at 66.9° < 0* < 70.5°. 
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* Experimental data (ref. 14) 

— Carey (eq. (5) and ref. 14) 

— Stephens (eq. (6) and ref. 19) 

— Present work (eq. (7)) 
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Figure 77. LIDCS plotted against transverse momentum for tt° production for E c 
6 * — 67° with data at 65.2° < 0* < 08.8°. 
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Figure 78. LIDCS plotted against transverse momentum for 7r° production for E cm — 30.8 GeV and 
53°. 
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Figure 79. FIDOS plotted against transverse momentum for 7T° production for E nil — 30.8 GeV and 
6* =. 90° with data at 30.6 < E rw <31 GeV. 



Pj_, GeV 

Figure 80. LIDOS plotted against transverse momentum for 7r° production for E cm = 45.1 GeV and 
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Figure 81. LIDCS plotted against transverse momentum for tt^ production for E cm = 45 GeV and#* = 90° 
with data at E rirt ~ 44.8 and 45.1 GeV. 
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Figure 82. LIDCS plotted against transverse momentum for 7T° production for E cm — 52.8 GeV and 
0* =r 11° with data at #* = 10° and 1 1°. 







Figure 83. LIDCS plotted against transverse momentum for 7r° production for E nr , — 52.8 GeV and 
0* = 16° with data at 0* = 15° and 17.5°. 



Figure 84. LIDCS plotted against transverse momentum for tt° production for E cm — 52.8 GeV and 
0* = 21° with data at 0* = 20° and 22°. 


65 






Figure 85. LIDCS plotted against transverse momentum for production for E nn = 52.8 GeV and 
6* = 7.5°. 



Figure 86. LIDCS plotted against transverse momentum for ir° production for E rm 
53°. 
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Figure 88. LIDCS plotted against transverse momentum for 7r° production for E cw — 7 GeV and 0 
with data at 58° < 0* < 61°. 
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Figure 90- LIDCS plotted against transverse momentum for tt 0 production for E cm — 62.6 GeV and 
$* = 90° with data at Zw, = 62.4 and 62.9 GeV. 
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Figure 91. LIDCS plotted against transverse momentum for tt° production for E rw = 9.8 GeV and 
0* ~ 45° with data at 44.9° < 0* < 45.6°. 



Figure 92. LIDCS plotted against transverse momentum for 7r° production for Erw 
0* = 64° with data at 63.6° < 0* < 65.1°. 
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Figure 93. LIDCS plotted against transverse momentum for tt° production for E rm — 9.8 GeV and 
0 * = 18° with data at 17.6° < 0* < 17.7°. 



Figure 94. LIDCS plotted against transverse momentum for 7 r° production for Ecru 
0* = 29° with data at 28.9° < 0* < 29.2°. 
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Figure 95. LIDCS plotted against transverse momentum for 7T° production for Ecm ~ 9.8 GeV and 
0 * - 55° with data at 54.6° < 0* < 55.7°. 



Figure 96. LIDCS plotted against transverse momentum for 7r° production for E rw 
0 * = 69° with data at 68.0° < 6>* < 69.9°. 
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Figure 97. LIDCS plotted against transverse momentum for tt~ production for Enn = 23 GeV and 
0 * = 90°. 
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Figure 98. LIDCS plotted against transverse momentum for 7 r“ production for E C m = 31 GeV and 
90°. 
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Figure 99. LIDCS plotted against transverse momentum for 7 r + production for E cm — 23 GeV and 

e* = 90°. 
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Figure 100. LIDCS plotted against transverse momentum for tt+ production for E rm = 31 GeV and 
90°. 
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Figure 102. LIDCS plotted against transverse momentum for i r + production for E cm — 45.0 GeV and 
0 * = 90° with data, at E cm = 45.0 and 44.8 GeV. 
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Figure 103. Parameter izations of total production cross section plotted with numerically integrated 
LTDCS parameterizations. 



Figure 104. Parameterizations of total tt production cross section plotted with numerically integrated 
LIDCS parameterizations. 
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Figure 105 . 7r° spectral distributions. For detailed comparison, t he horizont al axis is plot ted both linearly 
(upper figure) and logarithmically (lower figure). 
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Figure 106. 7T° spectral distributions. For comparison, horizontal axis is plotted both linearly (upper 
figure) and logarithmically (lower figure). 
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Figure 107. n spectral distributions. For comparison, horizontal axis is plotted both linearly (upper 
figure) and logarithmically (lower figure). 





Figure 108. 7r + spectral distribution. For 
and logarithmically (lower figure). 
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